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Abstract: We examine numerically different zero-dimensional reaction-diffusion processes 
as candidate toy models for high-energy QCD evolution. Of the models examined - Reggeon 
Field Theory, Directed Percolation and Reversible Processes ~ only the latter shows the 
behaviour commonly expected, namely an increase of the scattering amplitude with increas- 
ing rapidity. Further, we find that increasing recombination terms, quantum loops and the 
heuristic inclusion of a running of the couplings, generically slow down the evolution. 
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1. Introduction 

The seminal work of Munier and Peschanski [1] established a solid link between high energy 
QCD evolution and the physics of statistical mechanics systems. In that and ensuing 
works [2,3], the connection between the Balitsky-Kovchegov (BK) equation [4,5] and the 
Fisher-Kolmogorov-Petrovsky-Piscounov (FKPP) equation [6,7], used to describe the time 
evolution of certain quantities in some statistical mechanics systems, was exploited to give 
a general and elegant derivation of the energy dependence of the saturation scale and of 
geometric scaling [8,9]. 

The BK equation is the mean field limit of the more general set of evolution equations 
- the B-JIMWLK^ equations [4, 10-16] - describing the evolution of QCD scattering am- 
plitudes at high energy, i.e. the dynamics of the Colour Glass Condensate (CGC). Most 
studies of high energy QCD evolution have been performed using the BK equation. The 
reasons behind such preference are clear. Not only is the BK equation considerably more 
tractable than the B-JIMWLK set, but also the differences between the full B-JIMWLK 
evolution scheme and its truncated BK version are, perhaps surprisingly, small. The de- 
tailed numerical evaluation of these differences was performed in [17] by recasting the full 
evolution as a random walk in some functional space, once again highlighting the possible 
deep connection between QCD evolution and statistical mechanics. 

The B-JIMWLK evolution equations were derived under the explicit assumption of 
the scattering of a dilute projectile on a dense target and should, therefore, be applicable 

^Balitsky-Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov-Kovner. 
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only for that physical situation. Such a limitation renders the CGC framework - either in 
its simplified BK version or for the more complete B-JIMWLK set - insufficient to address 
the experimentally important collisions of two nuclei. Prom a more theoretical perspective, 
this incompleteness of B-JIMWLK has been shown to eventually lead to the violation of 
unitary constraints [18,19]. 

The inclusion of the missing ingredients - variously referred to as pomeron loops, 
fluctuations or wave function saturation effects - in an extended B-JIMWLK scheme has 
been attempted by several groups. A rather stringent constraint on the form of the complete 
evolution is imposed by what has been called the Dense-dilute Duality (DdD) [20,21]. This 
result, which establishes the equivalence of the dilute projectile, dense target case addressed 
by B-JIMWLK and its mirror image of a dense projectile and a dilute target, implies that 
the complete evolution kernel should necessarily be self-dual under the DdD. 

Of particular relevance to the subject matter of this paper are those attempts [22-35] 
to include the missing effects which resort to possible generalizations of the correspondence 
of high energy QCD and statistical mechanics beyond the proven instance of [1]. It has 
been argued [22] that the incompleteness of B-JIMWLK can be related to the breakdown 
of the FKPP description in systems involving a small number of objects. Thus, a powerful 
tool for the generalization of B-JIMWLK could be provided by insights from statistical 
mechanics systems where the discreteness effects that violate the FKPP description can be 
accounted for via the inclusion of stochastic effects. 

In the context of statistical mechanics systems, the failing FKPP description is reme- 
died by the inclusion of a stochastic (noise) term accounting for the effects of number 
fluctuations which become increasingly important with decreasing number of intervening 
objects. Different choices for the stochastic term to be added to the FKPP equation will, 
in general, lead to different dynamical behaviours. The choice of the noise term is thus 
a crucial step in obtaining a suitable description of high energy QCD evolution in the 
framework of stochastic processes. 

In general, there are two alternative strategies that can be followed to determine the 
form of the noise term. One can attempt to extend the known B-JIMWLK evolution 
equations by including the missing effects and then interpret the resulting equations in the 
stochastic language. This approach was followed in [25,28] where the obtained noise terms 
were, unfortunately, far too complex to allow any computation - analytical or numerical 
- to be carried out. Instead, one can choose to make an educated guess as to what the 
noise should be. Ultimately, for any such choice, the behaviour of the resulting Langevin 
equation will need to be compared with our general expectations for high energy evolution. 
In this paper we follow this latter strategy. 

Our analysis is performed in the general framework of reaction-diffusion processes, in 
particular in the Hamiltonian language of [36], which is equivalent to the formulation in 
terms of Langevin equations. A Langevin equation in which the stochastic term is such 
that the functional form of the two-point noise correlation function replicates that of the 
non-diffusive part of the FKPP equation is usually referred to as the stochastic FKPP 
(sFKPP) equation and describes the dynamics of a reaction-diffusion process belonging 
[37] to the universality class of the reversible processes. If the noise is chosen to have 
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a two-point correlator simply proportional to the quantity obeying the FKPP equation, 
the resulting Langevin equation belongs to a different universality class, that of directed 
percolation. Although other choices for the form of the stochastic term are certainly 
possible, we concentrate on these two since they have been previously considered as possible 
descriptions of high energy QCD in the framework of stochastic processes. 

We explore numerically the universality classes of reaction-diffusion processes with 
zero transverse dimensions. The Hamiltonian approach of [36] allows us to address both 
reversible processes and directed percolation universality classes within a single formula- 
tion. Further, the same formulation is well suited to enable a discussion of Reggeon Field 
Theory - both in its original form and with the extra quartic vertex recently considered 
in the context of high energy QCD evolution [34, 38] - which, although related to directed 
percolation [39], does not correspond to a well defined reaction-diffusion process. Finally, 
it also allows for a rather straightforward heuristic introduction of 'running coupling ef- 
fects'. Zero-dimensional toy models are thus used as a test bed for the far more involved 
realistic case with two transverse dimensions. By comparing the dynamical behaviour of 
the reaction-diffusion models with general expectations for high energy QCD evolution - 
geometrical scaling in BK, increase of the scattering amplitude with increasing rapidity, 
slowdown of evolution once couplings are allowed to run, diffusive scaling once stochas- 
tic effects are included - wc narrow down the number of possible candidate models. The 
central result of this paper is that, at least in relation to toy models with zero transverse 
dimensions, the only universality class of reaction-diffusion processes compatible with high 
energy QCD evolution is that of reversible processes. 

In Sec. 2 we introduce the formalism of [36], the toy models to be studied, an heuristic 
implementation of running of the couplings, and the classical solution of the problem. Our 
results are presented in Sec. 3 and a concluding discussion in Sec. 4. 

2. Formalism 

In this Section we start by presenting the basic formalism that will be used to discuss the 
different classes of reaction-diffusion models. Then we motivate a choice of parameters 
in the models, discuss an heuristic introduction of the running of the couplings, and the 
classical solutions to the problem. 

2.1 Basic setup 

We use the formalism developed in [36]. In this reference, the diffusive term in a stochastic 
equation plays the role of a kinetic term in a Hamiltonian formulation, while the details 
of the process are encoded in a reaction Hamiltonian which corresponds to a potential. 
Following the notation of [36], the reaction Hamiltonian for an elementary reaction kA A 
mA, where A is the strength of the reaction, reads 

HR{p,q) = ^{p^-p^)q\ (2.1) 

with [q,p] = 1. 
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Instead of discussing arbitrary cases, we focus in this paper in Hamiltonians with two 
powers of q at most. While higher vertices might be required by an eventual 'true' theory of 
high energy QCD, we restrict to this form which contains the vertices most often discussed 
in the literature. To pass to a more usual notation, we make the changes p —q, q ^ p, 
Hji —H, arriving at 

H{p, q) = aipq - a2qp^ - olz(^p + a4gV • (2-2) 

This expression corresponds to the Hamiltonian which rules, for reactions involving vertices 
up to fourth order, the evolution in rapidity y (the logarithm of the energy) of an auxiliary 
function F{y,q) [38], 

^E^ = -H{p,q)F{y,q). (2.3) 
dy 

The evolution starts from some functional form for the coupling to the projectile at y = 0, 
e.g. an eikonal coupling 

F(y = 0,g) = l-exp(-5ig), (2.4) 

where gi is the coupling parameter with the projectile. The relation with the transition 
amplitude ^/i(y) is given through 

iAfi{y)=F{y,q = gf), (2.5) 

where 5/ is the coupling parameter with the target. 

In (2.2) the variable p plays the role of a conjugate variable to q, p = —d/dq = —dq. 
It is, in the formalism of [36], the auxiliary degree of freedom required to study the non- 
equilibrium problem (for other approaches see [40-43]). Further, q{p) can be interpreted 
as the creation (annihilation) operator of the objects exchanged in the scattering. Within 
this interpretation, the property of the Hamiltonian (2.2), 

Hip, dp) = Hy^^y (^q, -—dq] , (2.6) 

amounts to a symmetry of the scattering process from the projectile and target points of 
view. This symmetry is strongly suggestive of the Dense-dilute Duality which has been 
postulated [20] to be an essential property of high energy QCD evolution Hamiltonians. 

In [36] the phase portrait determined by the zero energy trajectories, Hr{P, q) = 0, has 
been used to classify reaction-diffusion processes into universality classes. We will focus in 
the following cases: 

1. Reggeon Field Theory (RFT): 04 = 0, 02 = 03- It contains no quartic vertex and 
the splitting and recombination strengths are equal. It has no correspondence to a 
reaction-diffusion process, and it displays the known phenomenon that the amplitude 
vanishes with increasing rapidity which is usually interpreted as tunneling [44-46]. 
The effect of the quartic vertex has been examined recently [34] due to its possi- 
ble connection to high energy QCD evolution. The zero energy lines determine a 
triangular phase portrait given by 

p = 0, g = 0, q=^+p. (2.7) 
012 
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2. Directed Percolation (DP): it is the case for the stochastic process with allowed 
reactions 1 ^ 0, 1 2 and 2^1, leading to ai = /U — A, as = /i and 02 = 0:4 = a. 
It contains a quartic vertex with the same strength as the recombination vertex, and 
its phase portrait is determined by 

p = 0, q = 0, q= — . (2.8) 

"2(1 +p) 

Renormalization group arguments lead to the conclusion that the quartic vertex is 
irrelevant at dimensions greater than 4 (see e.g. [47]). On the other hand, the third 
trajectory in (2.8) goes into a straight line when expanded around p = 0. In both 
cases one returns to RFT, a fact that we will use below to fix the parameters for our 
numerical study. 

3. Reversible processes (RP): it is the case for the stochastic process with reactions 

1-^2 and 2^1, leading to ai = = fi and 02 = = cr. It contains a quartic 
vertex with the same strength as the recombination vertex and its phase portrait is 
determined by 

P=0, q = 0, P=-l, q=—. (2.9) 

This phase portrait cannot be deformed to a triangle and it does, therefore, belong to 
a different universality class than that of DP. The Langevin equation corresponding 

to this reaction-difussion problem is the stochastic Fisher-Kolmogorov-Pctrovsky-Pis- 
counov (see recent analyses in [48,49]), which has been proposed to play a role in 
high energy QCD evolution beyond JIMWLK [23,25,27,35]. 



2.2 Choice of parameters 

Instead of scanning a large region of the parameter space, we focus on the effect of the 
quartic vertex and of the difference between different reaction-diffusion processes. Thus, 
for RFT we choose values ai = 1 (this parameter can always be absorbed in a redefinition 
of rapidity) and 0:2 = 0:3 = 0.5. We have checked that with this set of parameters the 
conclusions we extract by numerically solving the evolution equation (2.3) in a restricted 
rapidity window < y < 5 hold for higher rapidities and are not a 'subasymptotic' effect, a 
situation which may happen for smaller values of 02 = as [38]. Besides, we will generically 
study the introduction of a quartic vertex in RFT as done in [34] though this situation, as 
RFT itself in zero dimensions, does not correspond to any reaction-diffusion process. 

For DP we adopt the following strategy: we approximate the phase portrait (2.8) 
by a triangle, which brings DP into RFT with a given quartic vertex, and identify the 
parameters with those in the phase portrait of RFT (2.7). This leads to the values ai = 1, 
a2 = 04 = 0.5 and 03 = 1.5. 

Finally, for RP we take ai = a3 = 1 and vary the remaining parameter in the region 
0.1 < 02 = 04 < 0.9. 
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2.3 Running of the couplings 

All couplings in these reaction-difussion processes are, at this level, fixed. Recently it has 
been proposed in the framework of a one-dimensional model, that effects of the running 
of the couplings may shift the contribution of effects of loops to higher rapidities [50] . To 
implement the running of the coupling in our approach we adopt an heuristic procedure. We 
interpret 1/q as some 'momentum' scale which determines a common logarithmic running 
for the four couplings in our Hamiltonian^ : 

<^i{q) = o^i'Pr^^ i = i,...,4, (2.10) 

for q < qo and frozen to the values aj discussed in the previous Subsection for q > qo- 
Q = lOqo plays the role of an inverse QCD scale, and we take two values of qo = 0.5 and 1. 
These choices are motivated by the need of making the effect of the running of the coupling 
noticeable in the restricted region of q we explore numerically. 

2.4 Classical solutions 

The Hamiltonian problem admits a classical solution: 

p = {-ai + a2p)p + 2{a3 - a4p)qp , (2.11) 
g = (ai - 2a2p)q + (-03 + 2a4p)g^ , 

with initial conditions 

q{y = 0)=9i, p{y = Y)=gf, (2.12) 

where Y is the total rapidity spanned between the projectile located at y = and the 
target located &t y = Y. Thus, the classical amplitude, i.e. the amplitude at tree level, 
can be computed standardly [46,52]: 

iAfriy) = 1 + 5; A.exp [-S{Y,qk,Pk)] , (2.13) 

k 

with the index k running over all possible solutions of the classical equations of motion 
(2.11)^, Ajfc = ±1 and S{Y,qk,pk) the action evaluated for the classical solutions. 



3. Results 



We begin this section with a description of the numerical method used to compute the 
transition amplitude. Following that, we study the rapidity evolution of the solutions 

^This choice of scale can be argued as follows: we identify the amplitude with that for the scattering 
of a dipole of size oc q, which is known to increase with increasing dipole size. Previous experience with 
running coupling effects in the BK equation [51] give us the hope that the concrete choice of scale will not 
alter the qualitative conclusions of this study. 

^Three solutions, one symmetric (A^ = +1) between projectile and target and two asymmetric {A^ = 
— 1), appear above some critical rapidity whose value depends on the chosen initial conditions, see a recent 
discussion in [52]. 
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for the fan case, RFT, DP and RP. Additionally, we examine the effects of an heuristic 
introduction of running of the couplings. Further, we phrase the effects of evolution in 
terms of a saturation scale. Finally, we compare the results of quantum and classical 
evolution. 



3.1 Numerical method 

To solve the differential equation (2.3) with evolution kernel (2.2) and initial conditions 
(2.4) which defines the quantum mechanical problem, we use a second order Rungc-Kutta 
method. We have discretized the gr-range in 500 points per unit which showed to be enough 
for a precision better than a few percent. The rapidity region studied is 

< y < 5 (3.1) 

(though we study values as high as 40 when comparing with the classical solution). The 
step in rapidity is correlated with the step in q [38]; we have used h = 6.25- 10~^. With this 
numerical method we obtain the values of Fi{y,q) in a g interval for different rapidities. 
The results depends also on the coupling constants: we use gi = 1 [38]. 

The solution to the classical equations of motion (2.11) with initial conditions (2.12) 
was obtained using the shooting method. After the classical trajectories were computed, 
the amplitude was calculated through (2.13) (the corresponding quantum expression is 
given by (2.5)). 

3.2 Numerical results 

We focus on the rapidity evolution of the solutions, 
(a) fan case 

First, as a check, we study the so-called fan case [53] in which there are no re- 
combination terms, from the projectile point of view, in the Hamiltonian (2.2), i.e. 
a2 = a4 = 0. The differential equation admits an analytical solution 



Ffan{y,q) = 1 - exp 



giqe^^y 



(3.2) 



We use this analytical result to perform a check of our numerical solution. As we 
can see in Fig. 1 there is a very good agreement between the analytical and the 
numerical results. Note that the BK equation sums the corresponding fan diagrams 
for BFKL pomerons [54] , leading a behaviour similar to what we found in this simple 
zero-dimensional case. 

(b) RFT 

Now we turn to RFT - Fig. 2 - which contains no quartic vertex. In the limit of 
high rapidity we find an exponential decay, a behaviour already found in [34, 38] and 
interpreted as a tunneling phenomenon [44-46]. 
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aj=l,a2=0,a3=0.5,a4=0 

I 1 1 1 1 1 1 

lines: numerical solution 



points: analytical solution 




Figure 1: Comparison of tlie analytical (symbols) and numerical results (lines) for the fan case, 
with gi ^ I. 

(c) DP 

The third case we study is DP, with the parameters fixed as discussed in Subsec- 
tion 2.2. What we find, Fig. 3, is that the evolution goes in the opposite direction 
to what is expected in high energy evolution, i.e. that the function moves to smaller 
values of q with increasing rapidity^. This poses serious doubts on DP as a candidate 
for a description of high energy QCD evolution. 

(d) RP 

Finally we turn to RP. As we see in Fig. 4 (left), the evolution leads the front 
towards smaller values of q as we increase the value of rapidity. So it behaves in the 
expected manner for high energy QCD evolution. The parameters 02,4, characterizing 
respectively the vertices 2 — > 1 and 2^2, are free. We find that increasing 02 = 04 
makes the evolution softer - Fig. 4 (right) -, again as expected for high energy 
evolution. 

(e) running coupling 

Now we turn to our heuristic implementation of the running of the couplings. In 
Fig. 5 we show the results for the fan case and for RP - the two cases where the 

^Our choice of parameters makes this behaviour visible in the y-range we study. Smaller values of 02,3,4 
make it noticeable only for larger values of y. 
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a^=l, a2=a3=0.5, a4=0 



0.8 



0.6 

LL 

0.4 



0.2 





0.5 1 1.5 2 

q 

Figure 2: Rapidity evolution of the solutions in the case of RFT. 

evolution goes in the expected direction - for the different freezing points for the 
couphng discussed in subsection 2.3. The evolution is clearly slowed down with the 
effect of the frozen procedure being more noticeable at the beginning of the evolution. 

(f) saturation scale 

To ensure that the effects of evolution in each case are more easily visible, we parallel 
the studies of the BK equation and introduce a 'saturation scale' l/qs{y), defined 
through the equation 

F{y,qs{y)) = K, (3.3) 

with K some fixed value of order 1/2. It is known [51] that the specific choice of the 
value of K leads to small changes in the leading y-behaviour of the saturation scale 
which, for the purposes of this study, are entirely negligible. In Fig. 6 we show the 
'saturation scale' obtained for the fan case, for DP and for RP. In the fan case, the 
running of the couplings slows down the increase of l/qs{y) (as in the BK equation). 
In DP I/qs (y) diminishes with increasing rapidity. Finally, in RP both the increase of 
the strengths of the vertices 2 — > 1 and 2 — > 2 and the running of the couplings slow 
down the increase of l/qs{y) - a behaviour which may be linked with the findings 
in [50] of a competition between the effects of recombination vertices, called 'Pomeron 
loops', and of running coupling. 

(g) classical solution 
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aj=l,a2=0.5,a3=1.5,a4=0.5 




Figure 3: Rapidity evolution of the solutions in the case of DP. 



ai=l, a2=0.3,a3=l, a4=0.3 ai=a3=l, y=5 




Figure 4: Rapidity evolution of the solutions in the case of RP for a2 = 0:4 = 0.3 (left). Solutions 
of RP for different values oi 02 ^ , at y — 5 (right). 

Finally, we want to compare the full quantum and the symmetrical classical solutions, 
see Subsection 2.4. This provides an evaluation of the effect of loops in the strict 
Quantum Field Theory sense. In Fig. 7 we show both the classical and quantum 
imaginary parts of the amplitude for RFT with a quartic vertex, for different values 
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Figure 5: Evolution of Fi{y,q) for given values of rapidity, comparing the behaviour with fixed 
(FC) and running (RC) couplings. Fan case (left). RP (right). 



of gi = Qf. The classical solution lies always above the quantum one, with the 
difference between both increasing with increasing rapidity. Similar conclusions are 
extracted in other cases. 

4. Conclusions 

In this paper we have examined numerically several reaction-diffusion processes in zero 
transverse dimensions as possible candidates for a toy model of high-energy QCD evolution. 
We have restricted ourselves to those models in which the Hamiltonian formulation [36] 
contain vertices with up to two derivatives: diffusion, splitting, recombination and 2^2, 
which allows a discussion of all common universality classes. We have also considered usual 
variants of such models which are not strictly reaction-diffusion processes. Their feasibility 
as toy models for high-energy QCD evolution has been evaluated in terms of the direction 
of the evolution with increasing evolution parameter, that is increasing rapidity. 

In the case of Reggeon Field Theory, which contains no quartic vertex and thus has 
no reaction-diffusion interpretation, we find the known behaviour of a vanishing amplitude 
with increasing rapidity. Directed Percolation leads to a behaviour of the solutions de- 
creasing with increasing rapidity. On the other hand, the limiting case of considering only 
fan diagrams (again with no reaction-diffusion counterpart) and Reversible Processes show 
solutions which increase with increasing rapidity. In the case of RP which, as DP, contains 
a quartic vertex, the evolution is slowed down by increasing the recombination and 2^2 
vertices. We have shown all these behaviours both at the level of the amplitudes and of an 
analogue of a saturation momentum. 

We have also introduced, in an heuristic manner, a running of the couplings, which 
tends to slow down the evolution in all cases, working in the same direction as the increase 
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aj=l,a2=0,a3=l,a4=0, Fj(y,q^)=0.5 ai=l,a2=0.5,a3=1.5,a4=0.5, Fi(y,q5)=0.2 




012345 012345 

y y 



aj=a3=l,Fj(y,q^)=0.5 



FC (a,=a4=0.1) ■ 

RC,<|).5(«2=«4=0-1) ■ 
RC,<,(a2=a4=0.1) 
FC(a,=a4=0.3) . 



Figure 6: l/qs{y) in the fan case for fixed and running couplings (top left). l/qs{y) in DP for 
fixed coupling (top right). \/qs{y) in RP with fixed coupling for two sets of a2 ^ on = 0.1 and 0.3, 
and for running coupling for a2 = 04 = 0.1 (bottom). 



of the recombination and 2 — > 2 vertices as found in [50]. Finally, we have numerically 
computed the classical solutions (those which sum tree diagrams), and found that they 
generically lie above the quantum ones, with their difference growing with increasing ra- 
pidity, so quantum effects tend to slow down the evolution. 

Summarizing, we have found that the only zero-dimensional reaction-diffusion process 
which shows a behaviour of the amplitude in agreement with the expectations from high 
energy QCD with rapidity is RP. This process is linked with sFKPP equation often used 
to discuss [23,25,27,35] high-energy QCD evolution. Increasing recombinations terms, and 
the inclusion of quantum loops and of a running of the coupling generically slow down the 
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ai=l, a2=a3=0.2, a4=0.04 



ai=l, a2=a3=0.2, a4=0.04 





ai=l, a2=a3=0.2, a4=0.04 




Figure 7: The classical symmetrical (dashed lines) and the quantum (solid line) imaginary part of 
the amplitude versus y. gi — gf ~ 0.1 (top left), = gf = 0.2 (top right), gi — gf = 1 (bottom). 



evolution. 
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